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The usefulness of glasses, and particularly of metallic glasses, in technological applications is 
often limited by their toughness, which is defined as the area under the stress vs. strain curve 
before plastic yielding. Recently toughness was found to increase significantly by the addition of 
small concentrations of foreign atoms that act as pinning centers. We model this phenomenon at 
zero temperature and quasi-static straining with randomly positioned particles that participate in 
the elastic deformation but are pinned in the non-afhne return to mechanical equilibrium. We find a 
very strong effect on toughness via the increase of both the shear modulus and the yield stress as a 
function of the density of pinned particles. Understanding the results calls for analyzing separately 
the elastic, or "Born term" and the contributions of the "excess modes" that result from glassy 
disorder. Finally we present a scaling theory that collapses the data on one universal curve as a 
function of rescaled variables. 



Introduction: In the field of metallic glasses an old 
but very important protocol called "micro-alloying" is 
getting recently very high attention. The protocol com- 
prises the addition of foreign atoms with concentration 
of the order of 1% to an otherwise regular metallic glass. 
Such micro-alloying results in dramatic changes in the 
physical and chemical properties of the material, includ- 
ing its glass forming ability diffusion properties 0, 
corrosion resistance [ij] and plastic properties These 
colossal changes, which are found even in the sub-percent 
alloying range as demonstrated by Ref. [1] for 0.5% Ag 
additions require a different approach compared to the 
normal alloying of two or more components in the mid- 
dle of the phase diagram. Changes in the electronic or 
chemical structure of the global composition are not suffi- 
cient to explain the measured large effects by such minor 
additions. This is also clearly shown in recent toughness 
measurements [H, where the addition of 1-2% of Si/Sn 
changes the notch toughness by 50%. 

In an earlier work Q it was suggested that the medium 
range order of the metallic glasses is reorganized and the 
stress distribution of that reconfiguration causes the 
change in long range property. In order to verify this 
suggestion we studied in a most recent MD-simulation 
the influence of pinned particles on the "global" prop- 
erty of the metallic alloy. Pinning particles in a very small 
concentration range act like too "big" or too "small" par- 
ticles that are introduced in a nearly perfect matrix of a 
glassy system. We found that the /3-process of the amor- 
phous system, which is well known to be responsible for 
aging and diffusion in the glassy state below the glass 
transition Q is totally suppressed by a small amount 
(up to 5%) of pinned particles. This is far below the per- 
colation limit of the system. In this Letter we show that 
the pinning protocol might provide a general approach to 
understand the effects of micro-alloying in bulk metallic 
glasses because it also affects the yield stress, the mod- 
ulus and the toughness. Predictions are made below for 



further comparison with experiments. 

Numerical simulations: To prepare accurate data 
for the present discussion we have performed 2- 
dimensional Molecular Dynamics simulations on a stan- 
dard binary system which is known to be a good glass 
former. We employed a 50 — 50 binary Lenard- Jones 
mixture with a potential energy for a pair of particles 
labeled i and j in the form 



a„(.,) = 4.,[(fJ)'^(fii)%A. + ,4,(I|) 

. (1) 



where the parameters Aq, Ai and A2 are added to 
smoothen the potential at a scaled cut-off oi r/a = 2.5 
upto the second derivative. The parameters cr^^, 
and cr^^ were chosen to be 2sin(7r/10), 2sin(7r/5) and 
1 respectively and e^^ = e^^ = 0.5, e^^ = l(see [^). 
The particle masses were taken to be equal. The sam- 
ples were prepared using high-temperature equilibration 
followed by a quench to zero temperature (T = 0.001) 
(see [l3l)- Our quench rate was 3.2 x 10^® in LJ units. 
All our procedures are performed with a fixed total num- 
ber of particles N and a constant volume such that the 
overall particles density p is strictly constant and equals 
0.976. 

Having prepared the amorphous solid in very much the 
standard way developed in the field, we now depart from 
the usual procedure by selecting randomly a set {P} of 
P particles that will be referred to below as the "pinned 
particles". These are chosen such that their density is 
Pp Pen where Pp and p„ are, respectively, the density 
of pinned and not pinned particles (p = p^ + p^). At this 
point we strain our sample using an athermal quasi-static 
(AQS) protocol to examine their stress vs. strain curves. 
In this procedure the particle positions in the system are 
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FIG. 1. (Color Online). Typical stress vs. strain curve ob- 
tained for AQS straining of one realization of a system of 
Af=20000 particles at different densities pp of pinned parti- 
cles. Note the significant change in the shear modulus (the 
slope at 7 = 0) and of the yield stress where the system yields 
to plastic flow. 



first changed by the affine transformation 



(2) 



During this affine step the pinned particles are allowed 
to participate like all the other particles. This trans- 
formation results in the system not being in mechanical 
equilibrium, and we therefore allow the non affine trans- 
formation r,; Ti + Ui which annuls the forces between 
the particles, returning the system to mechanical equi- 
librium. In this step the pinned particles are not allowed 
to participate; they are held fixed at the positions that 
they attained at the end of the affine step. The other 
particles are moving, directed by a gradient energy mini- 
mization, to seek new positions that annul the forces be- 
tween all the particles, pinned and not pinned. In other 
words during the energy minimization we add counter 
force such that the net force experienced by the pinned 
particles is constrained to be zero. Thus the main effect 
of the pinned particles is to constrain the system not to 
follow the standard conjugate gradient path but to find 
new "inherent states" or local minima in the energy land- 
scapes, different from those that could be found if all the 
particles were allowed to move in the non-affine step. 

In Fig. [l]we show the raw data of stress vs strain for 
this material with varying densities of pinned particles 
as presented in the inset. Shown are averages over 5 
realizations of the glass, with the same configuration of 
randomly chosen pinned particles. We observe that both 
the shear modulus and the yield stress (where the sys- 
tem yields to plastic flow) increase significantly when the 
density of pinned particles Pp is increased. In Fig. [2] we 
present the shear modulus as a function of Pp (see the 
middle curve with inverted triangles). In the rest of the 



FIG. 2. (Color Online). Shear modulus as a function of the 
density pp/po of pinned particles (N = 20000). In black in- 
verted triangles we show p as measured from the slope of Fig. 
[T] at 7 = 0. In red dots we show p as computed from the 
theory. The Born contribution is given by the black squares, 
showing an almost independence of pp. The non-affine con- 
tribution is shown in terms of the magenta circles. Of course 
the shear modulus itself is the difference between the Born 
term and the non-affine contribution. 



Letter we want to clarify this phenomenon in a quanti- 
tative way. 

Theory: By definition the shear modulus is the second 
derivative of the energy of the system with respect to the 
strain 7, i.e. 



1 (fU{ri, ■■■ ,rjv;7) 



V 



(3) 



In our process the full derivative with respect to 7 trans- 
lates to two contribution, one the direct partial derivative 
with respect to 7 and the other, via the chain rule, the 
contribution due to the non-affinc part of the transfor- 
mation: 



d d d dui _ d 

dl 97 .^7^ dui dj ~ 97 



y A.^ (4) 



where the second equality follows from the form of the 
non-affine transformation where dri = dui. To under- 
stand how to compute the constrained sum we need to 
recall that the force fi = —dU/dvi is zero before and 
after the affine and non-affine steps. Thus 
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Denote now as usual the Hessian matrix H and the non 



affine "force" H |17| : 
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Inverting Eq. ([5]) we find 



dui 



(7) 



Applying these results to the definition of /i we end up 
with the exact expression 



1 a^[/(ri,--. ,r^;7) 1 



V 



Qj2 



- ^ -i?^/ -^j , (8) 



where the first term is the well known Born contribution 
which we denote below as [Lb- The second term exists 
only due to the non-affine displacement Ui and it includes 
an unusual matrix of rank (iV — P) x N. Needless to say, 
before we compute the non-affinc contribution in Eq. (|5]) 
we need to remove the two Goldstone modes with A = 
which are the result of translation symmetry. 

It is very important to stress at this point that the sep- 
aration between the Born term and the non-affine term 
is not an arbitrary one. The Born term is very insensi- 
tive to pinned particles in our example, and this is usu- 
ally the case: it is only sensitive to average properties 
like density, average number of neighbors and interac- 
tions [l3|. In Fig. [5] we show the result of calculating 
the Born term for all our samples as a function of the 
density of pinned particles, (see black squares in Fig. [2), 
and there is only minor dependence. This is not the case 
for the non-affine term, whose direct calculation is also 
shown in the same figure in magenta circles. We see that 
this term changes significantly, taking upon itself the full 
blame of the change in the shear modulus as a function 
of the density of pinned particles. The difference of the 
two terms agrees to very high accuracy with the direct 
measurement of the shear modulus from the slopes of 
the curves in Fig. [T] at 7 = 0. We note that in other 
cases in which, say, the shear modulus depended on the 
quench rate from the melt, the physics was determined 
by the density of states which varied considerably with 
the quench rate [31 . This is NOT the case here. The 
full Hessian matrix is oblivious of the pinned particles, 
and its spectrum remains invariant when we change their 
density. The full reason for the change in shear modulus 
is embodied in the restricted sum in Eq. ([5]), respecting 
the fact that the pinned particles are not involved in the 
non-affine step of our straining protocol. 

Scaling theory: Having at hand the exact micro- 
scopic theory to evaluate the effect of the pinned particles 
on the shear modulus, we turn now to a scaling approach 
which may yield the phenomenological dependence of fi 
on Pp . We start with the scaling ansatz 



KPp) - m(Pp = 0) = e^^PpT{pp/po) 



(9) 



where the product ^abPp exhausts the dimensionality 
of the LHS and is a scaling function of a dimension- 
less variable. Expecting that the scaling function might 
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FIG. 3. (Color Online). Log-log plot of the rescaled shear 
modulus p{pp) — p{pp = K^abPp) ^ function of the 
dimensionless ratio pp/ po. The slope is o? = —0.25 ±0.02 and 
the data agree with the scaling function ansatz as represented 
in Eq. 



be represented by a power law we plot in a log-log plot 
p{pp) — p{pp = 0)/{£abPp) as a function of Pp/po- This 
plot is shown in Fig. [3l indicating that to a very good 
approximation the scaling function J- is indeed a power 
law. 



J^iPp/Po) fx (Pp/po)", a 



-0.25 ±0.02 



(10) 



collecting these results together indicates a phenomeno- 
logical law for the shear modulus in the form 



PiPp) - PiPp = 0) oc p^ 



0.75 



(11) 



It should be very interesting to compare this theoretical 
prediction to experiments on micro-alloying. 

Encouraged by this result we proceed further to seek a 
scaling form for the whole stress vs strain response as a 
function of rescaled variables with the aim of providing 
scaling collapse for the data in Fig. [TJ We start with the 
obvious scaling ansatz 



<^ii^Pp) ~ <^h.Pp = 0) 



3PpGil^Pp/Po) ■ (12) 



From the analysis of scaling of the shear modulus we 
conclude that 



lira Gil, Pp/po) oc p^ 

7— s-0 



-0.25 



(13) 



The exponent in the scaling law changes when 7 and a 
increase. We analyzed the data in the region of 7 « 7^ 
and found that 



lim g{-f,Pp/po) oc p^ 



-0.50 



(14) 



We can combine the last two equations into an interpola- 
tion formula which provides the global scaling form of Q. 
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exact theory matched perfectly with the simulation re- 
sults. Probably the most important result of this letter 
is the offered scaling theory for both the shear modulus 
and the entire stress vs. strain curves. This theory allows 
for a data collapse summarizing large amount of data in 
terms of a few simple functions of dimensional variables. 
It would be extremely interesting to test the applicability 
of similar scaling ideas to experimental examples of the 
effect of micro-alloying. 
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FIG. 4. (Color Online). The universal function if (7) which 
is obtained as explained in the text. This scaling function 
represents a collapse of all the data that appear in Fig. [1] 
Note that in this collapsed form the values of a^, as well as 
of become independent of the pinning. 



Defining cr^, as the value of a at which the stress reaches 
a maximum in Fig [TJ we write 



(15) 

where 1^(7) is expected to be a universal scaling function 
of 7, independent of concentration of pinned particles. 
Denoting 



^y{PpIpq)' 



a+[a^~a){pjp,)-^-^^ ' 



(16) 



we present in Fig. U] the LHS of Eq. ([T5|) divided by Q as 
a function of 7. The analysis indicates that the resulting 
plot should be the universal function ^^(7). 

Conclusions: We have demonstrated that our model 
micro-alloying which is realized by pinning randomly a 
small percentage of particles during the non-affine relax- 
ation step has a huge effect on the mechanical properties, 
i.e., shear modulus, yield-stress and toughness. For ex- 
ample the yield stress changes by more than 50% upon in- 
creasing the concentration of the pinned particles by 5%. 
The effect on the shear modulus had been explained by 
an exact theory separating the Born term from the non- 
affine contribution. It was shown that the latter takes 
full responsibility for the changes in shear modulus. The 
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